Entanglement detection in hybrid optomechanical systems 
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We study a device formed by a Bose-Einstein condensate (BEC) coupled to the field of a cavity with a 
moving end-mirror and find a working point such that the mirror-light entanglement is reproduced by the BEC- 
light quantum correlations. This provides an experimentally viable tool for inferring mirror-light entanglement 
with only a limited set of assumptions. We prove the existence of tripartite entanglement in the hybrid device, 
persisting up to temperatures of a few milli-Kelvin, and discuss a scheme to detect it. 
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I. INTRODUCTION 

Our ability to experimentally enforce quantum mechanical 
behaviors in complex systems of a diverse nature is improv- 
ing at a steady pace [ 1 ] and is swiftly leading us toward the 
grasping of the long craved quantumness at the macro scale. 
However, the faithful direct inference of quantum properties 
of parts of a complex systems or of a meso-system or macro- 
system is not a simple task to accomplish. This could be due 
to the intrinsic impossibility to directly address the system we 
are interested in or to its inherent complexity, which might 
face the lack of stringent and certifiable criteria for non clas- 
sicality. The idea of indirectly probing the system of interest 
by coupling it to fully controllable detection devices appears 
very appealing |2|,|3|]. It has spurred interest in designing tech- 
niques for the extraction of information from noisy or only 
partial data-sets gathered through observations of the detec- 
tion system. The main problem in such an approach is em- 
bodied by the set of extra assumptions that one has to make 
on the mechanisms to test and the retrodictive nature of the 
claims that can be made. In fact, one can a posteriori in- 
terpret a data-set and infer the properties of a system that is 
difficult to address by post-processing the outcomes of a few 
measurements. However, this requires assuming knowledge 
of the working principles of the effect that we want to probe. 

In this sense, it would be highly desirable to conceive in- 
direct detection strategies providing a faithful picture of the 
property to test by, for instance, "copying" it onto the de- 
tection stage, which would then be subjected to a direct es- 
timate process. By focussing on a hybrid optomechanical de- 
vice comprising an atomic, mechanical and optical mode, we 
propose a diagnostic tool for entanglement that operates along 
the lines of such desiderata. We show that the entanglement 
established between the mechanical and optical modes by the 
means of radiation-pressure coupling can indeed be "written" 
into the light-atom subsystem and directly read from it by 
means of standard, experimentally friendly homodyning. Our 
analysis goes further to reveal that such a possibility arises 
from the non-trivial entanglement sharing properties of the 
system at hand, which is indeed genuinely tripartite entangled. 
Very interestingly, we find that while any bipartite quantum 




FIG. 1: Sketch of the detection scheme. The probe field has the 
same periodicity of the cavity field. It impinges on a semi-reflecting 
mirror, creating a standing wave. Part of the probe light is transmitted 
and then measured with the output light of the cavity. 



correlation is bound to disappear at very low temperature, the 
multipartite content persists longer against the operating con- 
ditions, as a result of entanglement monogamy relations. Our 
analytical characterization considers all the relevant sources 
of detrimental effects in the system, is explicitly designed to 
be readily implemented in the lab and provides a first step into 
the assessment of multipartite entanglement in a macro-scale 
quantum system of enormous experimental interest. 



This paper is organized as follows. In section [TT] we de- 
scribe the physical system we consider and the equations of 
motion for the relevant variables. In section [III] we show the 
results for the stationary entanglement shared by the cavity 
field, the BEC and the optomechanical oscillator. In section 
IIVI we discuss how the entanglement between the three differ- 
ent subsystems can be experimentally measured and, finally, 
in section |V1 we summarize and conclude. 
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II. THE PHYSICAL MODEL 

The system we consider, sketched in Fig. Q] is made of 
a cavity where one of the two mirrors is free to move in 
the direction parallel to the cavity axis. In absence of radi- 
ation pressure from the cavity, the mirror moves harmonically 
with a frequency co m and is at equilibrium at a temperature 
T. The cavity is pumped through the fixed mirror by a laser 
of strength rj and is coupled to an intra-cavity BEC. We as- 
sume the latter in the weakly interacting regime U4J where we 
can expand the field operator in a classical part (representing 
the condensate wave-function) and a quantum part represent- 
ing quantum fluctuations. These are conveniently expressed 
in terms of Bogoliubov modes. It has been recently shown 
experimentally [5, 6] that the only Bogoliubov mode which 
is significantly coupled to the cavity is the one with momen- 
tum 2k c where k c is the cavity mode momentum. We also 
assume ultralow temperature of the condensate so that ther- 
mal fluctuations are negligible and the BEC wave-function is 
not affected by the coupling to the cavity field. We write the 
Hamiltonian of the system made of the cavity field, the mov- 
able mirror and the BEC as 



<H = tiM + ti c +<H A +<H A c + tiMc 



(1) 



where the mirror, cavity and BEC Hamiltonians are given re- 
spectively by 

= \ m ^ + ^ V 
Ha = hQ.c f c, (3) 
c Hq = S(wc — a>L)a^a — ihrj(a — a ' ). (4) 

Here, q and p are the mirror displacement and momentum op- 
erators while m and cj„, are its effective mass and oscillation 
frequency. The cavity (pump laser) frequency is u>c (u>l) and 
a is the cavity mode annihilation operator. The term propor- 
tional to 77= y/2Kfl/hciJi describes the pumping of the cavity 
and depends on the laser power % and the cavity decay rate 
k. Finally, Q and c are the frequency and bosonic annihilation 
operator of the Bogoliubov mode [cfr. Fig [T] . For small mir- 
ror displacements, the mirror-cavity interaction Hamiltonian 
can be written as I2B: 



(5) 



It is useful to introduce the quadrature operators for the Bo- 
goliubov mode 



G = 



c+c< 



P = i- 



V2 ' V2 ' 

so that the BEC-cavity interaction reads 

< Hac—^- — -a} a + h^Oa'a, 
2A a ^ 



(6) 



(7) 



where g is the vacuum Rabi frequency, iV is the condensate 
fraction, A a is the detuning of the atomic transition from the 



cavity frequency and the coupling rate f depends on the Bo- 
goliubov mode-function QSJ] - The same model can be realized 
with two BECs trapped in the same cavity, where the me- 
chanical oscillator is replaced by a second BEC. This setup 
can be realized straightforwardly as an extension of current 
experiments JUS- It is also worth stressing that, although 
our model shares, at first sight, some similarities with the 
one considered in Ref. II 1011 . it is distinctively different. First, 
the atomic media utilized in the two cases are different, with 
Ref. |H considering the bosonized version of the collective- 
spin operator of an atomic ensemble. Second, the resulting 
coupling Hamiltonians are built from quite distinct physical 
mechanisms. Third, the working conditions to be used in the 
two models in order to achieve interesting structures of shared 
quantum correlations are very different. We will comment on 
this latter point later. 

Any realistic description of the problem at hand should in- 
clude the most relevant sources of noise affecting the overall 
device, that is energy leakage from the cavity and mechani- 
cal Brownian motion at temperature T. The resulting open- 
system dynamics can be efficiently described by means of a 
set of Langevin equations for the vector of operators 



4>'=(x,y,q,P, Q,P), 



where 



a+a' 



V2 ' 



y = i- 



V2 



(8) 



are the position- and momentum-like quadrature operators of 
the cavity field. Under the experimentally reasonable assump- 
tion of an intense pumping field, we expand the relevant op- 
erators around the respective classical mean values (p s j as 
4'i=4's,j+54>j7 where 5^>/s are zero-mean fluctuation operators. 
The mean values are easily determined by solving the steady- 
state Langevin equations, which readily lead to 



Q s =- 



(9) 
(10) 



with 



kv = 



A 2 +K 2 



(that also allows us to determine x s and y s ) and the effective 
cavity detuning 



A-ojc-ojl 



g 2 No 
' 2A„ " 



(ii) 



In the rest of the paper, we will use A as a reference parame- 
ter, while all the other parameters entering the expression of A 
remain unspecified. At this point of our analysis, it is conve- 
nient to define dimensionless position and momentum opera- 
tors for the mirror as q- y/h/m(L> m q and p- ~\/Hma) m p and cast 
the Langevin equations for the fluctuations into the form 



d,54>=K5(t)+N, 



(12) 
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(a) (b) (c) 




FIG. 2: (a) Logarithmic negativity E A c against A (in units of <jj„,/27r=3xl0 6 s~') and the cavity-mirror coupling rate^. (b) Same as in panel 
(a) but for E M c- We have taken T=10 /uK, mechanical quality factor 3xl0 4 , m=50 ng, R=50 mW, cavity finesse F = 10 4 and (=x- The cavity 
decay rate is given by k = nc/lLF where c is the speed of light and the cavity length is L = 1 mm (c) Comparison between E AC (solid) and 
E M c (dotted) against T for A=2a) m and x=100 s . Other parameters are as in (a). 



where we have defined the vector of fluctuations 

64> T =(6x, 6y, 6q, dp, 6Q, 5P), 

the noise vector 

-\Jhmco m 

and the coupling matrix 



K = 



' -K 
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(13) 



that is readily determined following Ref. [8] and depends on 
the scaled coupling parameter \-x' ^/ti/mu) m and the mirror 
dissipation rate y. The noise vector contains the cavity input- 
noise operators 



OXifi — 



5ci;„+5a. 



Sym = i- 



V2 
5a)—5ah 



V2 



(14) 
(15) 



such that {5ai n )-0 and {5ai„(t)5a in (t'))=5{t-t'). Moreover, we 
have introduced the zero-mean Langevin force operator | re- 
sponsible for Brownian motion of the mechanical system at 
temperature T. Small values of y and low operating tempera- 
tures allow to retain an effective Markovian description of the 
mechanical Brownian motion, where ifTTll 



<^)£(/')>-^<5(f-0 

PB 



(16) 



with PB-fili^kBT) and kg being the Boltzmann constant. 
Moreover, we have checked the validity of the first-order ex- 
pansion of the operators (which holds only for small quantum 



fluctuations) by looking at the stability of the steady-state so- 
lution to the Langevin equations, which is guaranteed when 7C 
has all negative eigenvalues. Using A as a free parameter, the 
solution is stable when A>0 and (x , < ^ K ^ which implies that 
the mirror-BEC-induced back-action is small compared to the 
cavity dissipation. Such conditions are met across our work. 



III. STATIONARY ENTANGLEMENT 

We can now address the entanglement-inference scheme 
that is the focus of our study. We focus on the f— >oo regime 
where any transient in the dynamics of the system has faded 
away and steady-state quantum correlations are set. The main 
tool of our assessment is the covariance matrix of the state of 
our hybrid system, which is defined as IU2I1 



(17) 



and contains full information on a Gaussian state such as ours. 
Using the Langevin equations one can show that *V fulfills the 
equation 



where 



and 



1C'V+'V<K T = -£>, 



£>=diag[*, k, 0, y(2h+ 1 ), 0, 0] 



l 



(18) 



exp(2/3 B w m )-l 



(19) 



is the thermal mean-occupation number of the mechanical 
state. This allows us to find *V for any values of the relevant 
parameters. We call A, M, C the atomic mode, the mirror one 
and cavity field respectively and quantify the entanglement 
between any two modes a and /3 (a,/3=A, M, C) using the log- 
arithmic negativity E aB =- log2v 11411 . Here, v is the small- 
est symplectic eigenvalue of the the matrix < V' a „-P r V a ^P , 
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P=diag(l, 1,1,-1) performs momentum-inversion in phase- 
space and "Yob is the reduced covariance matrix of a and 

In the following two subsections, we discuss separately the 
bipartite and the multipartite cases. 



A. Bipartite entanglement 

Here we analyze the stationary entanglement in the three 
possible bipartitions of the system. Let us start with the sym- 
metric situation where Q.-u) m and where the cavity-mirror 
coupling equals the cavity-atoms coupling ((=x)- We always 
restrict ourselves to a stable regime for the linear Langevin 
system. In this case and for very low initial temperature of 
the mechanical mode (we take r=10yuK), we find that entan- 
glement is generated in the stationary state between the Bo- 
goliubov mode and the cavity field as well within the field- 
mechanical mode subsystem. Due to the symmetric coupling 
Eac is very similar to Emc, as shown in Fig. |2a-b). This 
provides an interesting diagnostic tool: the inaccessibility of 
the mirror mode makes the inference of the optomechanical 
entanglement a hard task needing cleverly designed, although 
experimentally challenging, indirect methods OLljl]. How- 
ever, in light of the recent demonstration of controllability of 
intra-cavity atomic systems JUGii, we can think of inferring 
the pure optomechanical entanglement simply by measuring 
the more accessible correlations between the optical field and 
the Bogoliubov mode. In the range of parameters considered 
here, the mirror-atoms entanglement Eam is always zero. This 
is in contrast with the results in Ref. Ill Oil , where the emer- 
gence of Eam is due to the use of an effective negative detun- 
ing that regulates the free evolution of the collective atomic 
quadrature. In our case, such evolution is ruled by the fre- 
quency Q, which is always positive. Therefore, the two mod- 
els can access quite different regions of the parameter-space. 
In turn, this implies that, while our system does not allow for 
hybrid atom-mirror entanglement, Genes et aim Ref. IToll did 
not achieve the symmetric Emc -Eac situation revealed above. 

We now study how Eac an d Emc decay with T. The re- 
sults plotted in Fig. [2] (c) show that the two entanglements 
are indistinguishable and therefore the BEC can still be used 
as an entanglement probe. As expected the entanglement de- 
cays when increasing the environment temperature and disap- 
pears for r>0.1mK. Higher critical temperatures for the dis- 
appearance of entanglement are found for larger cavity quality 
factors and more intense pumps. We consider two different 
regimes, where we relax the symmetric conditions between 
the mirror and the Bogoliubov mode. In the first situation, 
we take co m =Cl and vary the coupling rates. In the second, 
we take symmetric couplings £=x and change the frequencies. 
For u> m =Cl, Eac an d Emc are strongly affected by the change 
in g and^f. As shown in Fig. [3] (a), Eac grows continuously 
for larger f while Emc decreases slightly. In this regime, it 
is clear that small inaccuracies of the ratio £/% do not affect 
the mirror-cavity entanglement, therefore confirming the role 
of the BEC as a minimal disturbance probe. A more involved 
situation occurs when we keep a>„, fixed and change £1 We 
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FIG. 3: (a) Logarithmic negativity E AC (solid) and E MC (dashed) 
against £/x for A=2co m , Other parameters are as in Fig. [2] (b) We 
show Eac (line) and E M c (dashed) against £l/u> m . The inset shows 
the entanglement functions around Cl=a> m for T=1/jK and F=4x 10 4 . 



take T-ljjK, cavity finesse F=4xl0 4 and find a sharp peak 
at Q.-(jj m where Emc-Eac, see [Fig. [3](b)]. While Emc in- 
creases slowly with Cl/a) m (the Bogoliubov mode goes out of 
resonance from the cavity and decouples from the rest of the 
system), E A c reaches its maximum at £2/w m =s2. These results 
show that by changing the Bogoliubov frequency, for exam- 
ple varying the longitudinal trapping frequency, the BEC acts 
as a switch for the mirror-cavity entanglement, inhibiting or 
enhancing it. 



B. Genuine multipartite entanglement 

We now study the existence of genuine multipartite entan- 
glement in our system looking at the logarithmic negativity in 
each one-vs-two-mode bipartition. According to Ref. II 1711 . if 
in a tripartite state all such bipartitions are inseparable, gen- 
uine multipartite entanglement is shared. For the systems at 
hand, we indeed find inseparability of the bipartitions A\MC, 
M\AC and C\AM [see Fig. 0), which strongly confirms the 
presence of tripartite entanglement up to very high temper- 
atures. Evidently, multipartite entanglement persists up to 
r~0.01K, which is much larger than the critical one for the 
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idea is to shine the BEC with a probe laser in a standing wave 
|[l8ll forming a small angle with respect to the cavity axis, as 
sketched in Fig. Q] For an off-resonant probe field with po- 
larization perpendicular to that of the primary cavity field, the 
additional Hamiltonian terms read 



< H P =% 



Ap+Up 



f 



d 3 rcos 2 (k c x+k z z)^ 



a p a P , 



(21) 



where Ap is the detuning of the probe field from the pump- 
ing laser, ap is the corresponding annihilation operator, 
Up-gpl Aap is the light-atom coupling constant, k z is the 
wave-vector of the probing field along a direction orthogonal 
to the cavity axis and 4* is the condensate field operator. For 
a BEC strongly confined in the plane orthogonal to the cav- 
ity axis, we can neglect excitations of transverse modes and 
factorize the field operator as 



FIG. 4: We show (i, j, k=A, M, C) and the genuine tripartite en- 
tanglement G, ri against T for Q=<±) m and;^=f . We have E A \ MC =E M \ AC 
(solid), Eqam (dashed) and G,„ (dot-dashed). Other parameters are 
as in Fig. [2] 



disappearance of bipartite entanglement (~8 x 10 5 K). Such 
a result spurs the quantitative study of the genuine tripartite- 
entanglement content of the overall state II 1511 . For pure mul- 
tipartite states, this is based on monogamy inequalities valid 
for the squared logarithmic negativity, which turns out to be a 
proper entanglement monotone. For mixed states, the convex- 
roof extension of such a measure is required, albeit restricted 
to the class of Gaussian states. More explicitly, we aim at 
determining 



V=Mx)f(y)fQd, 



(22) 



Gtn= min[G,| ;it -G, 7 -G,i] 



(20) 



with TL(ijk) the permutation of indices i, j,k=A, M,C and G 
the convex roof of the squared logarithmic negativity for a 
given bipartition of the system. In general, the evaluation of 
Gtn is very demanding. However, for o) m =Q andx=(, the co- 
variance matrix *V is symmetric under the permutation of A 
and M, which greatly simplifies the calculations: The residual 
tripartite entanglement can be thus determined against the ef- 
fects of T, showing a non monotonic behavior. To understand 
this, we take i=C, j-M and k-A (any other combination could 
be considered). Similar to E^c and Emc, Gc\a,m decays very 
quickly as T increases, thus biasing the competition between 
Gqma and G C |a + Gqm- Analogously to Eqam, G C \a + G C \m 
decreases very slowly with T, thus determining an overall in- 
crease of the residual entanglement. However, as T is raised 
further, the system tends toward two-mode biseparability and 
any tripartite entanglement is washed out [see Fig. [4). 



IV. PROBING EXPERIMENTALLY THE STATIONARY 
ENTANGLEMENT 

We can now discuss the direct observation of the entangle- 
ment between the Bogoliubov mode and the cavity field. The 



where f(y) is a function peaked around the cavity axis and 
characterized by a width cr. We then expand ifr(x) as [4] 

iff(x)= yfN^i// (x)+ [ u ko-(x)c ka --vl a .(x)cl r \. (23) 

jfc>0,cr=± 

We assume that the width of the BEC in the transverse di- 
rection is much smaller than the periodicity of the probe field 
along the z direction. This implies that transverse excitations 
are suppressed. Neglecting terms 0(k 2 cr 2 ) we get 

( Hap-%[a p +— ^^+( P Ci\alap-iHT]p(ap-& f p ), (24) 

where £p is the same as £ with Uo—*Up. The Langevin equa- 
tion for the probe field reads 



d t ap=T]p-i 



a N oUp , A 
Ap+——+£ P Q 



dp+ y2i<api n —Kap 



(25) 



Assuming again an intense pump field, the steady intensity of 
the probe field is 



\ap\ 



Al+K 2 ' 



(26) 



where 



A a ^0 r n 

Ap=Ap+——+^ P Q s 



We also assume that £><kQ, £ and ap<^as ■ The equation of 
motion for the fluctuations of the probe field is 

d t 6ap— — iAp5cip—i£pa p 5Q+ yf2K~Sapj n —Kp5ap. (27) 

In order to map the Bogoliuobov mode onto the probe field, 
and following the same technique as in 111311 . we choose 
Ap=£l»Kp,£pap. The equations of motion for the slowly 
varying variables d{t)-e' a ' d(t) thus read 



(pap ~ / — ~ 
o t 5a- - i — —c+ \2K5api„-Kp5ap. 
V2 



(28) 
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If the decay of the probe cavity is faster than the dynamics of 
the Bogoliubov modes, the probe field follows adiabatically 
the dynamics of the latter. Using the cavity input-output rela- 
tions 1H2I1 we get 

Sa Pout =-i^^-c+6api„ (29) 

which shows how the Bogoliubov mode is mapped onto the 
output probe field. To measure the entanglement between the 
field of the primary cavity and the Bogoliubov mode one can 
homodyne the cavity field and rotate the quadrature of the 
probe. Analogously, one could map q onto a further field so 
that the light-matter correlations are changed into amenable 
light-light ones 11311 . 

V. CONCLUSIONS 

We have demonstrated the proof of principle of the use of a 
BEC as a diagnostic tool to determine the elusive mirror-light 



entanglement in a hybrid optomechanical device. No extra 
assumption is required for such an inference protocol that can 
be demonstrated using current experimental abilities JHl 0]. 
The entanglement structure within the system is intriguing and 
includes the case of tripartite inseparability, resilient to tem- 
peratures of up to a few milli-Kelvin and fully experimentally 
characterizable using existing technology. This demonstrates 
the viability of our scheme for diagnostics which contributes 
to the study of a complex three-body interaction exhibiting 
genuine quantum features at the macroscale. 
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